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A search has been performed for photons originating in the decay of a neutral long-lived particle, 
exploiting the capabilities of the ATLAS electromagnetic calorimeter to make precise measurements 
of the flight direction of photons, as well as the calorimeter's excellent time resolution. The search has 
been made in the diphoton plus missing transverse energy final state, using the full data sample of 
4.8 fb _1 of 7 TeV proton-proton collisions collected in 201 1 with the ATLAS detector at the LHC. No 
excess is observed above the background expected from Standard Model processes. The results are 
used to set exclusion limits in the context of Gauge Mediated Supersymmetry Breaking models, with 
the lightest neutralino being the next-to-lightest supersymmetric particle and decaying with a lifetime 
in excess of 0.25 ns into a photon and a gravitino. 
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I. INTRODUCTION 

Supersymmetry (SUSY) [lHslj a theoretically well- 
motivated candidate for physics beyond the Standard 
Model (SM), predicts the existence of a new SUSY part- 
ner (sparticle) for each of the SM particles, with iden- 
tical quantum numbers except differing by half a unit 
of spin. In R-parity conserving SUSY models |10l - ll4 |] 
, these sparticles could be produced in pairs in proton- 
proton (pp) collisions at the CERN Large Hadron Col- 
lider (LHC), and would decay in cascades involving other 
sparticles and SM particles until the lightest SUSY par- 
ticle (LSP), which is stable, is produced. In Gauge Medi- 
ated Supersymmetry Breaking (GMSB) models [15l - l20| . 
the LSP is the gravitino (67). GMSB phenomenology 
is largely determined by the properties of the next-to- 
lightest supersymmetric particle (NLSP). The results of 
this analysis are presented in the context of the so-called 
Snowmass Points and Slopes parameter set 8 (SPS8) |2l| . 
which describes a set of minimal GMSB models with the 
lightest neutralino (Xi) as the NLSP. In the SPS8 set of 
models, the effective scale of SUSY breaking, denoted A, 
is a free parameter. In addition, the Xi proper decay 
length, ct(Xi), is a free parameter of the theory. 

For the range of SPS8 A values relevant for this anal- 
ysis, SUSY production is dominated by electroweak pair 
production of gauginos, and in particular of X2X1 and 
Xi Xi pairs. In the SPS8 models, the dominant de- 
cay mode of the NLSP is Xi ->■ 7 + G. Previous AT- 
LAS analyses have assumed prompt NLSP decays, with 
cr(Xi) < 0.1 mm, and therefore searched for an excess 
production of diphoton events that, due to the escaping 
gravitinos, exhibit significant missing transverse momen- 
tum (Kj? iss ). The latest such ATLAS results [H use 
the full 2011 dataset and, within the context of SPS8 
models, exclude values of A < 196 TeV, corresponding to 
m(X?) > 280 GeV, at the 95% confidence level (CL). The 
limits on SPS8 models are less stringent in the case of a 
longer-lived NLSP. For example, recent CMS 95% CL 
limits 23] , obtained using the E™ ms spectrum of events 



with at least three jets and one or two photons, coupled 
with measurements of the photon arrival time, require 
m(X?) > 220 GeV for cr(Xi) values up to 500 mm. 

The analysis reported in this paper uses the full data 
sample of 4.8 fb _1 of 7 TeV pp collisions collected in 2011 
with the ATLAS detector at the LHC, and considers the 
scenario where the Xi has a finite lifetime and can travel 
some distance from its production point before decay- 
ing. The search is performed in an inclusive sample of 
candidate diphoton + E™ lss events. 

The long-lived NLSP scenario introduces the possibil- 
ity of a decay photon being produced after a finite de- 
lay and with a flight direction that does not point back 
to the primary vertex (PV) of the event. The analy- 
sis searches for such "non-pointing photons" by exploit- 
ing the fine segmentation of the ATLAS electromagnetic 
(EM) calorimeter to measure the flight direction of pho- 
tons. The variable used as a measure of the degree of 
non-pointing of the photon is zdca, the difference be- 
tween the z-coordinatqj of the photon extrapolated back 
to its distance-of-closest-approach (DCA) to the beam- 
line (i.e. x — y — 0) and zpy, the z-coordinate of the PV. 
The search for non-pointing photons is then performed 
by fitting the shape of the zdca distribution obtained for 
photons in the signal region, defined as diphoton events 
with E™ lss > 75 GeV, to a combination of templates that 
describe the zdca distribution for the expected signal 
and background events. In addition, the excellent time 
resolution of the calorimeter is exploited to measure the 
arrival times of the photons, providing a cross-check of 
the results. 



1 ATLAS uses a right-handed coordinate system with its origin 
at the nominal interaction point (IP) in the center of the de- 
tector and the z-axis along the beam pipe. The i-axis points 
from the IP to the center of the LHC ring, and the j/-axis points 
upward. Cylindrical coordinates (r, tfi) are used in the trans- 
verse plane, <f> being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 9 as 
r\ = -lntan(0/2). 



II. THE ATLAS DETECTOR 

The ATLAS detector [2J] covers nearly the entire solid 
angle around the collision point, and consists of an in- 
ner tracking detector surrounded by a solenoid, EM and 
hadronic calorimeters, and a muon spectrometer incor- 
porating three large toroidal magnet systems. The AT- 
LAS inner-detector system (ID) is immersed in a 2 T 
axial magnetic field, provided by a thin superconducting 
solenoid located before the calorimeters, and provides 
charged particle tracking in the pseudorapidity range 
| r\ | < 2.5. The ID consists of three detector subsys- 
tems, beginning closest to the beamline with the high- 
granularity silicon pixel detector, followed at larger radii 
by the silicon microstrip tracker and then the straw-tube- 
based transition radiation tracker. The ID allows an ac- 
curate reconstruction of tracks from the primary pp colli- 
sion and precise determination of the location of the PV. 
The ID also identifies tracks from secondary vertices, per- 
mitting the efficient identification of photons that convert 
to electron-positron pairs as they pass through the de- 
tector material. 

This analysis relies heavily on the capabilities of the 
ATLAS calorimeter system, which covers the pseudora- 
pidity range \r]\ < 4.9. Finely-segmented EM calorime- 
try is provided by barrel (|?7| < 1.475) and end-cap 
(1.375 < M < 3.2) lead/liquid-argon (lead/LAr) EM 
sampling calorimeters. An additional thin LAr pre- 
sampler covering |?/| < 1.8 allows corrections for en- 
ergy losses in material upstream of the EM calorimeters. 
Hadronic calorimctry is provided by a steel/scintillator- 
tile calorimeter, segmented into three barrel structures 
within |?7| < 1.7, and two copper/LAr hadronic end- 
cap calorimeters. The solid angle coverage is completed 
with forward copper/LAr and tungsten/LAr calorimeter 
modules optimized for EM and hadronic measurements, 
respectively. Outside the ATLAS calorimeters lies the 
muon spectrometer, which identifies and measures the 
deflection of muons up to 1 77 1 = 2.7, in a magnetic field 
generated by superconducting air-core toroidal magnet 
systems. 



A. Pointing resolution 

The EM calorimeter is segmented into three layers in 
depth, that are used to measure the longitudinal profile of 
the shower. The first layer uses highly granular "strips" 
segmented in the rj direction, designed to allow efficient 
discrimination between single photon showers and two 
overlapping showers originating, for example, from the 
decay of a n° meson. The second layer collects most 
of the energy deposited in the calorimeter by EM show- 
ers initiated by electrons or photons. Very high energy 
showers can leave significant energy deposits in the third 
layer, which can also be used to correct for energy leak- 
age beyond the EM calorimeter. By measuring precisely 
the centroids of the EM shower in the first and second 



layers, the flight direction of photons can be determined. 
In the ATLAS H — > 77 analysis [25[ that contributed 
to the discovery of a Higgs-like particle, this capability 
of the EM calorimeter was used to help choose the PV 
from which the two photons originated, thereby improv- 
ing the diphoton invariant mass resolution and sensitivity 
of the search. The analysis described in this paper uses 
the measurement of the photon flight direction to search 
for photons that do not point back to the PV. The angu- 
lar resolution of the EM calorimeter's measurement of the 
flight direction of prompt photons is about 60 mr&d/y/E, 
where E is the energy measured in GeV. This angular 
precision corresponds, in the EM barrel calorimeter, to 
a resolution on zdca of about 15 mm for prompt pho- 
tons with energies in the range of 50-100 GeV. Given 
the geometry, the zdca resolution is worse for photons 
reconstructed in the end-cap calorimeters, so the point- 
ing analysis is restricted to photon candidates in the EM 
barrel calorimeter. 

While the geometry of the EM calorimeter has been 
optimized for detecting particles that point back to near 
the nominal interaction point at the center of the detector 
(i.e. x = y = z = 0), the fine segmentation allows good 
pointing performance to be achieved over a wide range of 
photon impact angles. Figure [1] shows, as a function of 
I £dca I, the expected pointing resolution for SPS8 signal 
photons from Monte Carlo (MC) simulations, obtained 
by fitting to a Gaussian the difference between the value 
of .zdca obtained from the calorimeter measurement and 
the MC generator-level information. The pointing resolu- 
tion degrades with increasing |zdca|, but remains much 
smaller than |zdca| in the region where the signal is ex- 
pected. 

The calorimeter pointing performance has been veri- 
fied in data by using the finite spread of the LHC colli- 
sion region along the z-axis. Superimposed on Figure[T]is 
the pointing resolution achieved for a sample of electrons 
from Z — > ee events, where the distance, zpv, between 
the PV and the nominal center of the detector serves 
the role of zdca- In this case, the pointing resolution is 
obtained by fitting to a Gaussian the difference between 
zpv, as determined with high precision using tracking in- 
formation, and the calorimeter measurement of the origin 
of the electron, along the beamline. Figure Q] shows that 
similar pointing performance is observed for photons and 
for electrons, as expected given their similar EM shower 
developments. This similarity validates the use of a sam- 
ple of electrons from Z — > ee events to study the pointing 
performance for photons. The expected pointing perfor- 
mance for electrons in a MC sample of Z — > ee events is 
also shown on Figure [T] and is consistent with the data. 
The level of agreement between MC simulation and data 
over the range of values that can be accessed in the data 
gives confidence in the extrapolation using MC simula- 
tion to the larger deviations from pointing characteristic 
of the signal photons. 
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FIG. 1. Pointing resolution obtained for EM showers in the ATLAS LAr EM barrel calorimeter. The pointing resolution for 
photons from GMSB signal MC samples is plotted as a function of |.zdca|. The pointing resolution is also shown for Z — > ee 
data and MC, for which the primary vertex position, 2pv, serves the role of 2dca- 



B. Timing resolution 

Photons from NLSP decays would reach the LAr 
calorimeter with a slight delay compared to prompt pho- 
tons. This delay results mostly from the flight time of the 
heavy NLSP, as well as some effect due to the longer ge- 
ometric path of a non-pointing photon produced in the 
NLSP decay. The EM calorimeter, with its novel "ac- 
cordion" design, and its readout, which incorporates fast 
shaping, has excellent timing performance. Quality con- 
trol tests during production of the electronics required 
the clock jitter on the LAr readout boards to be less 
than 20 ps, with typical values of 10 ps [26j. Calibra- 
tion tests of the overall electronic readout performed in 
situ in the ATLAS cavern show a timing resolution of 
« 70 ps [23], limited not by the readout but by the jit- 
ter of the calibration pulse injection system. Test-beam 
measurements 28] of production EM barrel calorimeter 
modules demonstrated a timing resolution of ss 100 ps in 
response to high energy electrons. 

For this analysis, the arrival time of an EM shower is 
measured using the second-layer EM calorimeter cell with 
the maximum energy deposit. During 2011, the various 
LAr channels were timed-in online with a precision of 
about 1 ns. A large sample of W — > ev events was used 



to determine several calibration corrections which were 
applied to further optimize the timing resolution for EM 
clusters. The calibration includes corrections of various 
offsets in the timing of individual channels, corrections 
for the energy dependence of the timing, and flight-path 
corrections depending on the position of the PV. The 
corrections determined using the W — > ev events were 
subsequently applied to electron candidates in Z — » ee 
events to validate the procedure as well as to determine 
the timing performance in an independent data sample. 
Figure [5] shows the time resolution achieved as a func- 
tion of the energy deposited in the second-layer cell used 
in the time measurement. The time resolution, a{t), is 
expected to follow the form u(t) — a/E b, where E is 
the energy measured in GeV and © indicates addition in 
quadrature. Superimposed on Figure [2] is the result of a 
fit of the resolution to this form, where the parameters 
a and b multiply the so-called noise term and constant 
term, respectively. A timing resolution of ~ 290 ps is 
achieved for a large energy deposit. By comparing the 
arrival times of the two electrons in Z — > ee events, this 
resolution is understood to include a correlated contri- 
bution of sw 220 ps, as expected due to the spread in pp 
collision times caused by the lengths of individual proton 
bunches along the LHC beamline. Subtracting this beam 



contribution in quadrature, the obtained timing resolu- 
tion for the calorimeter is « 190 ps. 
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FIG. 2. Time resolution obtained for EM showers in the AT- 
LAS LAr EM barrel calorimeter, as a function of the en- 
ergy deposited in the second-layer cell with the maximum 
deposited energy. Superimposed is the result of the fit de- 
scribed in the text. The data are shown for electrons read 
out using high gain, and the errors shown are statistical only. 

To cover the full dynamic range of physics signals of 
interest, the ATLAS LAr calorimeter readout boards [26| 
employ three overlapping linear gain scales, dubbed high, 
medium, and low, where the relative gain of successive 
scale is reduced by a factor of about ten. For a given 
event, any individual LAr readout channel is digitized 
using the gain scale that provides optimal energy resolu- 
tion, given the energy deposited in that calorimeter cell. 
The results in Figure [2] are those obtained for electrons 
where the time was measured using a second-layer cell 
read out using high gain, for which the W — > ev sample 
used to calibrate the timing is large. Calibration samples 
for the medium and low gain scales are smaller, resulting 
in reduced precision. The timing resolutions obtained 
are about 400 ps for medium gain and about 1 ns for low 
gain. 



While all background studies, apart from some cross- 
checks, are performed with data, MC simulations are 
used to study the response to SPS8 GMSB signal mod- 
els, as a function of the free parameters A and ct. The 
other SPS8 model parameters are fixed to the following 
values: the messenger mass M mcss = 2A, the number of 
SU(5) messengers N 5 — 1, the ratio of the vacuum expec- 
tation values of the two Higgs doublets tan/3 = 15, and 
the Higgs-sector mixing parameter \x > [21j . The SPS8 
SUSY mass spectra, branching ratios and decay widths 
are calculated using ISAJET [29[ version 7.80. The sig- 
nal yield is normalized to the central value of the GMSB 
signal cross section, as calculated at next-to-leading or- 
der (NLO) using PROSPINO [30| version 2.1 with the 
CTEQ6.6m [H| parton distribution functions (PDFs). 
The total uncertainty on the signal cross section, times 
the signal acceptance and efficiency, from varying the 
PDFs and factorization and renormalization scales, as de- 
scribed in Ref. [32J . ranges from 4.7% to 6.4%, depending 
on A. 

The Herwig-|--|- generator version 2.4.2 [33] with 
MRST2007LO* [H PDFs is used to generate the signal 
MC samples. The branching ratio of the X\ — > 7 + G de- 
cay mode is fixed to unity. Signal samples were generated 
for fixed values of A, ranging from 70 TeV to 210 TeV in 
10 TeV steps, and for a variety of values of the NLSP life- 
time. The performance for any lifetime can be obtained 
by reweighting as a function of lifetime the existing sam- 
ples for a given value of A. The current analysis considers 
NLSP lifetime values in excess of 0.25 ns. 

All MC samples were processed with the simula- 
tion [13 of the ATLAS detector based on GEANT4 [H 
and were reconstructed with the same algorithms used 
for the data. The presence of additional pp interactions 
(pile-up) as a function of the instantaneous luminosity 
is taken into account by overlaying simulated minimum 
bias events according to the distribution of the number 
of pile-up interactions observed in data, with an aver- 
age of about nine interactions per bunch crossing. The 
MC samples were generated with a spatial distribution 
of the PV consistent with that observed in the data. The 
simulated LAr time values are smeared such that the 
time resolution matches the performance, as described in 
Sec. Ill Bl observed in data. 



IV. OBJECT RECONSTRUCTION AND 
IDENTIFICATION 



III. DATA AND MONTE CARLO SIMULATION 

SAMPLES 

This analysis uses a data sample of pp collisions at a 
center-of-mass-energy of ^/s = 7 TeV, recorded with the 
ATLAS detector in 2011. The data sample, after apply- 
ing quality criteria that require all ATLAS subdetector 
systems to be functioning normally, corresponds to a to- 
tal integrated luminosity of 4.8 ±0.1 fb -1 . 



The reconstruction of converted and unconverted pho- 
tons and of electrons is described in Refs. [37| and [38| . 
respectively. Shape variables computed from the lateral 
and longitudinal energy profiles of the EM shower in the 
calorimeter are used to identify photons and discrimi- 
nate against backgrounds. Two sets of selection crite- 
ria, denoted "loose" and "tight", are defined [37]. The 
loose photon identification, designed for high photon effi- 
ciency with modest background rejection, uses variables 



describing the shower shape in the second layer of the EM 
calorimeter, as well as leakage of energy into the hadronic 
calorimeter. The tight photon identification, designed for 
higher purity photon identification with still reasonable 
efficiency, includes more stringent cuts on these variables, 
as well as requirements on additional variables describ- 
ing the shower in the hrst layer of the EM calorimeter. 
The various selection criteria do not depend on the trans- 
verse momentum of the photon (Et), but do vary as a 
function of r\ in order to take into account variations in 
the calorimeter geometry and in the thickness of the up- 
stream material. For more details, see Ref. [37j. 

The measurement of E!f lss [39| is based on energy de- 
posits in calorimeter cells inside three-dimensional clus- 
ters with 1 77 1 < 4.5. The energy of each cluster is cali- 
brated to correct for the different response to electromag- 
netically and hadronically induced showers, energy loss 
in dead material, and out-of-cluster energy. The value 
of E™ lss is corrected for contributions from any muons 
identified in the event, by adding in the energy derived 
from the properties of reconstructed muon tracks. 



V. EVENT SELECTION 

The selected events were collected by an online trigger 
requiring the presence of at least two loose photon candi- 
dates, each with £t > 20 GeV and \r]\ < 2.5. To ensure 
the selected events resulted from a beam collision, events 
are required to have at least one PV candidate with five 
or more associated tracks. In case of multiple vertices, 
the PV is chosen as the vertex with the greatest sum of 
the square of the transverse momenta of all associated 
tracks. 

The offline photon selection requires the two photon 
candidates to each have 3? > 50 GeV, and to satisfy 
\tj\ < 2.37, excluding the transition region of 1.37 < \tj\ < 
1.52 between the barrel and end-cap EM calorimeters. 
In addition, both photons are required to be isolated: 
after correcting for contributions from pile-up and the 
deposition ascribed to the photon itself [37|, the trans- 
verse energy deposited in the calorimeter in a cone of 
AR = ^/(A?/) 2 + (A0) 2 =0.2 around each photon can- 
didate must be less than 5 GeV [40{. 

Due mostly to the inclusion of cuts on the EM shower 
shape in the very finely segmented strips of the first EM 
calorimeter layer, the efficiency of the tight photon iden- 
tification decreases with |zdca| for values of |zdca| larger 
than « 100 mm. The loose efficiency remains flat over a 
wider range of |zdca|, up to values of « 250 mm, after 
which it decreases less rapidly than the tight efficiency. 
The event selection requires at least one of the isolated 
photon candidates to pass the tight photon identification 
requirements, while the other must pass the loose photon 
identification cuts. The selected sample will therefore be 
referred to hereafter as the Tight-Loose (TL) diphoton 
sample. To reduce the potential bias in the pointing mea- 
surement that results from applying the photon identifi- 



cation requirements, only the loose photon in each event 
is examined for evidence of non-pointing. The 77 cut on 
the loose photon is tightened to restrict it to lie in the 
EM barrel calorimeter, namely \rj\ < 1.37. 

The TL diphoton sample is divided into exclusive sub- 
samples according to the value of E!f lss . The TL sample 
with E^ iBS < 20 GeV is used, as described in Sec. IVTBl 
to model the prompt backgrounds. The TL events with 
intermediate E™ iss values, namely 20 GeV < E!f iss < 
75 GeV, are used as a control sample to validate the 
analysis procedure. The final signal region, which con- 
tains a total of 46 selected events, is defined by applying 
to the TL diphoton sample the additional requirement 
that the value of E™ [ss exceeds 75 GeV. 

Table U summarizes the total acceptance times effi- 
ciency of the selection requirements, for examples of SPS8 
signal model points with various values of A and r. For 
fixed A, the acceptance falls approximately exponentially 
with increasing r, dominated by the requirement that 
both NLSPs decay inside the ATLAS tracking detector 
(which extends to a radius of 107 cm) so that the de- 
cay photons are detected by the EM calorimeters. For 
fixed t, the acceptance increases with increasing A, since 
the SUSY particle masses increase, leading the decay 
cascades to produce, on average, higher E!f lss and also 
higher E^ values of the decay photons. 



T 




A (TeV) 




(ns) 


80 


120 


160 


0.25 


15.3 ±0.3 


29.6 ±0.3 


45.1 ±0.3 


1 


11.1 ±0.1 


27.0 ±0.2 


35.9 ±0.3 


6 


2.01 ±0.02 


5.38 ± 0.02 


8.06 ± 0.06 


20 


0.39 ±0.01 


1.006 ± 0.005 


1.43 ±0.01 


40 


0.175 ±0.005 


0.384 ± 0.002 


0.510 ±0.004 


80 


0.090 ± 0.004 


0.164 ±0.001 


0.196 ±0.002 



TABLE I. The total signal acceptance times efficiency, in per- 
cent, of the event selection requirements, for sample SPS8 
model points with various A and r values. The uncertainties 
shown are statistical only. 



VI. SIGNAL AND BACKGROUND MODELING 



A. SPS8 GMSB signal 

The shape of the zdca distribution, alternately de- 
noted hereafter as the pointing distribution, that is ex- 
pected for photons from NLSP decays in events passing 
the selection cuts is determined using the SPS8 GMSB 
MC signal samples described previously, for various val- 
ues of A and r. The signal pointing distributions, nor- 
malized to unit area, are used as signal templates, here- 
after referred to as T s ; g . Since the T s j g shape is deter- 
mined using MC simulations, systematic uncertainties 



in the shape are included to account for possible differ- 
ences in pointing performance between MC simulation 
and data. In particular, the presence of pile-up in the col- 
lisions could impact the pointing resolution, due both to 
energy deposits in the calorimeter from additional mini- 
mum bias collisions and to the possibility of misidentify- 
ing the correct PV. These effects are modeled in the MC 
simulation. However, as a conservative estimate of the 
systematic T s - lg shape variations that could occur due to 
pile-up, the shapes from the two subsamples are obtained 
by dividing the MC samples roughly in two according to 
the number of PV candidates identified in each event. 
The T S i g distributions, along with their statistical and to- 
tal uncertainties, are shown in Figure [3] for A = 120 TeV 
and for NLSP lifetime values of r = 0.5 and 30 ns. 



B. Backgrounds 

The background is expected to be completely domi- 
nated by pp collision events, with possible backgrounds 
due to cosmic rays, beam-halo events, or other non- 
collision processes being negligible. The source of the 
loose photon in background events contributing to the 
selected TL sample is expected to be either a prompt 
photon, an electron misidentified as a photon, or a jet 
misidentihed as a photon. In each case, the object pro- 
viding the loose photon signature originates from the PV. 
The pointing and timing distributions expected for these 
background sources are determined using data control 
samples. 

Given their similar EM shower developments, the 
pointing resolution is similar for prompt photons and for 
electrons. The zdca distribution expected for prompt 
photons and for electrons is therefore modeled using elec- 
trons in Z — > ee data events. The Z — > ee event selection 
requires a pair of oppositely charged electron candidates, 
each of which has px > 25 GeV and I77I < 2.37 (exclud- 
ing the transition region between the barrel and end-cap 
calorimeters). One of the electron candidates, dubbed 
the "tag", is required to pass additional selection cuts 
on the EM shower topology and tracking information 
designed to identify a high-purity sample of electrons. 
The second electron candidate, dubbed the "probe", is 
restricted to the range |?7| < 1.37. The dielectron in- 
variant mass is required to agree with the value of the Z 
mass within 10 GeV, a requirement that produces a suffi- 
ciently clean sample of Z — > ee events. To avoid any bias, 
the pointing resolution for electrons is determined using 
the distribution measured for the probe electrons. The 
pointing distribution determined from Z — > ee events, 
normalized to unit area, is used as the pointing template 
for prompt photons and electrons, and is referred to here- 
after as T e /j. 

While the EM showers of electrons and photons 
are similar, there are some differences. In particu- 
lar, electrons traversing the material of the ID may 
emit bremsstrahlung photons, widening the resulting EM 



shower. In addition, photons can convert into electron- 
positron pairs in the material of the ID. In general, the 
EM showers of unconverted photons are slightly narrower 
than those of electrons, which are in turn slightly nar- 
rower than those of converted photons. The EM com- 
ponent of the background in the signal region includes a 
mixture of electrons, converted photons, and unconverted 
photons. Therefore, using electrons from Z — > ee events 
to model the EM showers of the loose photon candidates 
in the signal region can slightly underestimate the point- 
ing resolution in some cases, and slightly overestimate 
it in others. The pointing distribution from Z — > ee 
events is taken as the nominal T e i^ shape and the ex- 
pected distributions from MC samples of unconverted 
and converted photons are separately taken to provide 
conservative estimates of the possible variations in the 
T e / 7 shape which could result from not separating these 
various contributions. The T e / 7 distribution, along with 
its statistical and total uncertainties, is shown superim- 
posed on Figure [31 

Due to their wider showers in the calorimeter, jets have 
a wider zdca distribution than prompt photons and elec- 
trons. The sample of events passing the TL selection, but 
with the additional requirement that E T 11SS < 20 GeV, is 
used as a data control sample that includes jets with 
properties similar to the background contributions ex- 
pected in the signal region. The E™ lss requirement serves 
to render negligible any possible signal contribution in 
this control sample. The shape of the zdca distribu- 
tion for the loose photon in these events, normalized to 
unit area, is used as a template, referred to hereafter as 
T E miss <2 o GeV' m the final fit to the signal region. The 
TL sample with E™ iss < 20 GeV should be dominated by 
jet-jet, jet-7 and 77 events. Therefore, the T E mis S<2 o GeV 
template includes contributions from photons as well as 
from jets faking the loose photon signature. When using 
the template in the fit to extract the final results, it is not 
necessary to separate the photon and jet contributions. 
Instead, the relative fraction of the two background tem- 
plates is treated as a nuisance parameter in the fitting 
procedure, as discussed in Sec. IVIII1 

The pointing resolution depends on the value of Et 
of the photon candidate. Applying the shape of the 
T B mis 8<20 GoV template to describe events in the signal 
region, defined with E^ lss > 75 GeV, therefore implic- 
itly relies on the assumption that the Et distributions for 
photon candidates are similar in both regions. However, 
since E™ 135 is essentially a negative vector sum of the Et 
values of the energy depositions in the calorimeter, it is 
expected that there should be a correlation between the 
value of E™ 185 and the Et distributions of the physics 
objects in the event. This correlation is indeed observed 
in the TL control regions: the samples with higher S™ lss 
values have higher average photon Et values. Increas- 
ing the minimum Et requirement on the photons to 
60 GeV in the E™ iss < 20 GeV control sample selects 
events with kinematic properties which are more simi- 
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FIG. 3. The zoca templates from Z — > ee events, from the TL control sample with _E™ 1SS less than 20 GeV, and for MC 
simulations of GMSB signals with A = 120 TeV and values for the NLSP lifetime of r = 0.5 and 30 ns. The data points show 
the statistical errors, while the shaded bands show the total uncertainties, with statistical and systematic contributions added 
in quadrature. The first (last) bin includes the contribution from underflows (overflows). 



lar to those of events in the signal region, and therefore 
the Et > 60 GeV cut is applied to determine the nomi- 



nal T K 



is <20 GeV 



template shape. Conservative possible 



systematic variations of this shape are taken to be the 
template shapes with Ef cuts of 50 GeV and 70 GeV 



on the photon candidates. The T E 



3 <20 GeV 



template, 



along with its statistical and total uncertainties, is shown 
superimposed on Figure [3] 

The distribution of arrival times of the photon candi- 
dates is used in the analysis as a cross-check. Figure 0] 
shows the timing distribution expected for selected signal 
events, for A = 120 TeV and for NLSP lifetime values of 
t = 0.5 and 30 ns. The expectations for the backgrounds 
are determined using the same data control samples de- 
scribed above. It is expected that the performance of 
the calorimeter timing measurement, as determined us- 
ing the second-layer cell with the maximum deposited 
energy, should be rather insensitive to the details of the 
EM shower development. It was verified that the tim- 
ing distribution of electrons in Z — > ee events is very 
similar to that of loose photon candidates in the TL con- 



trol sample with E™ 1SS < 20 GeV. Therefore, the timing 
distribution determined with the larger Z — > ee sample 
is characteristic of the timing performance expected for 
all prompt backgrounds, and is shown superimposed on 
Figure |U 



VII. SYSTEMATIC UNCERTAINTIES 

When fitting the photon pointing distribution in data 
to the templates describing the expectations from signal 
and background, the total number of events is normalized 
to the 46 events observed in the data in the signal region. 
In addition, the background templates are determined 
using data. Thus, the analysis does not rely on predic- 
tions of the background normalization or composition, 
and there are therefore no systematic uncertainties to 
consider regarding the normalization of the backgrounds. 
As a result, the various systematic uncertainties relevant 
for the analysis can be divided into two types, namely 
"flat" uncertainties that are not a function of zdca but 
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FIG. 4. The distribution of photon arrival times (t 7 ) expected 
for SPS8 GMSB signal models with A = 120 TeV and for 
NLSP lifetime values of r = 0.5 and 30 ns. Superimposed is 
the expectation for prompt backgrounds, as determined using 
electrons from Z — > ee events. The uncertainties shown are 
statistical only. 



determination of the diphoton trigger efficiency using a 
bootstrap method [42|, and a contribution of ±2% taken 
as a conservative estimate of any possible dependence 
of the trigger efficiency on zdca- Uncertainties on the 
photon selection, the photon energy scale, and the de- 
tailed material composition of the detector, as described 
in Ref. [22|], result in an uncertainty of ±4.4%. The un- 
certainty due to the photon isolation requirement was 
estimated by varying the energy leakage and the pile-up 
corrections independently, resulting in an uncertainty of 
±1.4%. Systematic uncertainties due to the E™ lss recon- 
struction, estimated by varying the cluster energies and 
the E™ lss resolution between the measured performance 
and MC expectations [39|, contribute an uncertainty in 
the range of ±(1.1-8.2)%, with the higher uncertainty 
values applicable for lower values of A. As described in 
Sec. lIIIl variations in the calculated NLO signal cross sec- 
tions, times the signal acceptance and efficiency, at the 
level of ±(4.7-6.4)% occur when varying the PDFs and 
factorization and renormalization scales. Adding these 
flat systematic uncertainties in quadrature gives a total 
systematic uncertainty on the signal yield in the range 
of ±(7.2-11.7)%, to which is added a contribution in the 
range of ±(0.7-5.0)% due to statistical uncertainties in 
the signal MC predictions. 



VIII. RESULTS AND INTERPRETATION 



that affect the overall expected signal yield, and "shape" 
uncertainties related to the shapes of the unit-normalized 
signal and background pointing templates. The shape 
uncertainties, which are correlated across zdca bins, are 
discussed in Sec. IVI[ their sizes within each bin are de- 
picted in Fig. [3] The various flat systematic uncertainties 
are summarized in Table UH and discussed in more detail 
below. 



Source of Uncertainty 


Value 


Integrated Luminosity 


±1.8% 


Trigger Efficiency 


±2.1% 


Photon ID and Et Scale/Resolution 


±4.4% 


Photon Isolation 


±1.4% 


E T " BB : E T Scale/Resolution 


±(1.1-8.2)% 


Signal PDF and Scale Uncertainties 


±(4.7-6.4)% 


Total Systematic Uncertainty 


±(7.2-11.7)% 



TABLE II. Summary of systematic uncertainties on the total 
signal yield. The final row provides the total of these system- 
atic uncertainties, calculated as the quadrature sum of the 
various contributions. There is an additional contribution, 
not shown in the Table, due to MC statistics, that ranges 
between ±0.7% and ±5.0%, depending mostly on the NLSP 
lifetime. 



The uncertainty 

±1.8% 



on the integrated luminosity is 
41] . The uncertainty on the trigger efficiency 



is ±2.1%, which includes a contribution of ±0.5% for the 



The zdca distribution for the 46 loose photons of the 
events in the signal region with _E™ 1SS > 75 GeV is shown 
in Fig. \5\ As expected for SM backgrounds, the distri- 
bution is rather narrow, and there is no obvious sign 
of a significant excess in the tails that would be ex- 
pected for GMSB signal photons originating in the de- 
cays of long-lived neutralinos. There are three events 
with |zdca| > 200 mm, including one with a value of 
zdca = ±752 mm. Some additional information about 
these three events is summarized in Table [TTT1 

The timing distribution for the 46 events in the signal 
region with E™ 1 ^ > 75 GeV is shown in Fig. [51 The tim- 
ing distribution is rather narrow, in agreement with the 
background-only expectation that is shown superimposed 
on Fig. [51 and there is no significant excess in the positive 
tail that would be expected for GMSB signal photons. 
The photon with the largest time value has t ~ 1.2 ns, 
and is measured using a calorimeter cell that was read 
out with medium gain. This photon corresponds to one 
of the three events with | zdca | > 200 mm that are de- 
scribed in Table [TTTl but not to the most extreme pointing 
outlier. For the other two events in Table [TTTT the timing 
is consistent with the hypothesis that the photon candi- 
date is in-time. 

To determine the final results, the pointing distribution 
shown in Fig. [5] is fitted using the pointing templates de- 
scribed previously. The binning shown in Fig. [5] is used 
in the fit. The background contribution is modeled in 
the fit as a weighted sum of the T e / 7 and Tg;miss <2 Q GcV 
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FIG. 5. The zdca distribution for the 46 loose photon candidates of the events in the signal region. Superimposed are the 
results of the background-only fit, as well as the results of the signal-plus-background fit for the case of A = 120 TeV and 
r = 6 ns. The hatching shows the total uncertainties in each bin for the signal-plus-background fit, for which the fitted signal 
strength is \x — 0.20 ± 0.19. The first (last) bin includes the contribution from underflows (overflows). The inlay shows an 
expanded view of the central region, near zdca — 0. 



Run 


Event 
Number 


E^ iBB Loose Photon 




Tight Photon 




Number 


(GeV) E T (GcV) zdca (mm) 


ty (ps) 


Et (GcV) zdca (mm) 


tj (ps) 


186721 


30399675 


77.1 75.9 -274 


360 


72.0 22 


580 


187552 


14929851 


77.3 59.4 -262 


1200 


87.2 -120 


240 


191920 


14157929 


77.9 56.6 752 


2 


54.2 5 


-200 



TABLE III. Some parameters of the three "outlier" events mentioned in the text. 



templates. Background-only fits are performed to deter- 
mine the compatibility of the observed pointing distribu- 
tion with the background-only hypothesis. Signal-plus- 
background fits are performed to determine, via profile 
likelihood fits, the 95% CL limit on the signal strength, 
(i, defined as the number of fitted signal events divided by 
the SPS8 expectation for the signal yield. In both cases, 
the overall normalization is constrained to the 46 events 
observed in the signal region in data, and the relative 
weighting of the two background templates is treated as 
a nuisance parameter in the fit. 



Table IIVI shows the number of observed events in 
the various bins of zdca shown in Fig. O except that 
the number of zdca bins is reduced for display pur- 
poses by a factor of two by taking the absolute value 
of zdca- Included in Table ITVl are the results from the 
background-only and signal-plus-background fit, for the 
case of A = 120 TeV and r = 6 ns. Comparing the 
first two rows of Table QV] shows that there is fair agree- 
ment between the observed data and the results of the 
background-only fit, though there is some excess seen 
in the data for larger |zdca| values. The po value for 



Fit 
Type 



Event 
Type" 



Range of |zdca| values [mm] 



0-20 20-40 40-60 60-80 



D-100 100-200 200-400 400-600 > 600 



Data 



27 



7 



1 







1 



10 



Bkg Only Bkg 25.0±2.2 9.1±0.8 3.8±0.3 2.1+0.5 1.4±0.4 3.0TT.1 1.3±0.5 0.2±0.1 0.08±0.03 

Signal Total 25.1±4.2 9.3±1.5 3.3±0.7 1.6±0.6 1.1±0.4 2.6±1.0 1.8±0.8 0.7±0.5 0.5±0.4 

Plus Sig 0.7±0.6 0.5±0.5 0.4±0.3 0.3±0.3 0.3±0.3 1.2±1.1 1.3±1.2 0.6±0.5 0.4±0.4 

Bkg Bkg 24.4±4.2 8.8±1.5 2.9T0.8 1.3+0.7 0.8±0.6 1.4TT.5 0.5±0.7 O.lzfcO.l 0.03±0.04 

TABLE IV. Integrals over various zdca ranges of the distributions shown in Fig. [5] for the 46 loose photon candidates in 
the signal region. The numbers of events observed in data are shown, as well as the results of a background-only fit and a 
signal-plus-background fit for the case of A = 120 TeV and r = 6 ns. The fitted signal strength is /i = 0.20 ± 0.19. The errors 
shown correspond to the sum of statistical and systematic uncertainties. The numbers of signal and background events from 
the signal-plus-background fit are negatively correlated. 
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FIG. 6. The distribution of arrival times (i 7 ) for the 46 loose 
photon candidates of the events in the signal region. Superim- 
posed for comparison is the shape of the timing distribution 
expected for background only, normalized to 46 total events. 



the background-only hypothesis is w 0.060, indicating 
that the slight excess has a significance equivalent to 
rj 1.5c. The excess is dominated by the outlier pho- 
ton with Zdca = 752 mm; removing this event from the 
distribution and performing a new fit, one obtains a pq 
value of ss 0.30, indicating in this case a much better 
agreement with the background-only model. 

The last three rows of Table HVl show the results of the 
signal-plus-background fit for the case of A = 120 TeV 
and t — 6 ns, including the total number of fitted events 
in each |zdca| bin, as well as the separate contributions 
from signal and from background. For this case, the 
signal-plus-background fit returns a central value for the 
signal strength of /i = 0.20±0.19. The fit results are used 



to determine, via the CLs method [43|, 95% CL limits on 
the number of signal events. The results, derived within 
the RooStats framework [44j , are determined for both the 
observed limit, where the fit is performed to the pointing 
distribution of the 46 observed events in the signal region, 
and for the expected limit, where the fit is performed to 
ensembles of pseudo-experiments generated according to 
the background-only hypothesis. The slight excess for 
larger |zdca| values leads to a result where the observed 
limit on the number of signal events is somewhat less 
restrictive than the expected limit. For the example of 
A = 120 TeV and r = 6 ns, the observed (expected) 95% 
CL limit is 18.3 (9.8) signal events. 

By repeating the statistical procedure for various A 
and t values, the limits are determined as a function of 
these SPS8 model parameters. The left plot of Fig. [7] 
shows the 95% CL limits on the number of signal events 
versus r, for the case with A = 120 TeV. The right plot 
of Fig. [7] shows the 95% CL limits on the allowed cross 
section versus r, also for the case with A = 120 TeV. 
Each plot includes a curve indicating the SPS8 theory 
prediction for A = 120 TeV. The intersections where the 
limits cross the theory prediction show that, for A = 
120 TeV, values of r below 8.7 ns are excluded at 95% 
CL, whereas the expected limit would exclude values of 
t below 14.6 ns. 

Comparing with the theoretical cross section of the 
SPS8 GMSB model, the results are converted into an ex- 
clusion region in the two-dimensional plane of r versus 
A, as shown in Fig. [5] Also shown in the figure are corre- 
sponding limits on the lifetime versus the masses of the 
lightest neutralino and lightest charginos, where the re- 
lation between A and sparticle masses is taken from the 
theory. 



IX. CONCLUSIONS 

A search has been performed for non-pointing photons 
in the diphoton plus E™ lss final state, using the full data 
sample of 7 TeV pp collisions recorded in 2011 with the 
ATLAS detector at the LHC. The analysis uses the ca- 
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FIG. 7. 95% CL limits on (left) the number of signal events and (right) the SPS8 signal cross section, as a function of NLSP 
lifetime, for the case of A = 120 TeV. The region below the limit curve is excluded at 95% CL. 
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FIG. 8. The expected and observed limits in the plane of NLSP lifetime versus A (and also versus the Xi or Xi masses), for 
the SPS8 model. Linear interpolations are shown to connect between A values, separated by 10 TeV, for which MC signal 
samples are available. The region excluded at 95% CL is shown as the blue hatched area. The limit is not shown below an 
NLSP lifetime of 0.25 ns, which, due to the MC signal samples available, is the smallest value considered in the analysis. 
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pability of the ATLAS LAr calorimeter to measure the 
flight direction of photons to perform the search, with 
the precision measurement of the arrival time of photons 
used as a cross-check of the results. 

No significant evidence for non-pointing photons is ob- 
served. Interpreted in the context of the GMSB SPS8 
benchmark model, the results provide 95% CL exclu- 
sion limits in the plane of r (the lifetime of the lightest 
neutralino) versus A (the effective scale of SUSY break- 
ing) or, alternatively, versus the mass of the lightest neu- 
tralino. For example, for A = 70 TeV (160 TeV), NLSP 
lifetimes between 0.25 and 50.7 ns (2.7 ns) are excluded 
at 95% CL. 
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